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Nanoelectronic  Modeling  Software  Development 
Contract  Number:  DAAH01-97-C-R276 


PROJECT  OBJECTIVE 

This  Phase  I  project  has  two  primary  objectives.  The  first 
objective  is  to  develop  an  engineering  model  of  Resonant 
Tunneling  Devices  (RTD)  to  be  implemented  in  AIM-Spice,  a 
commercial  software  package  running  under  Windows  95/NT.  The 
development  of  a  reliable  device  model  is  a  prerequisite  for 
the  design  and  production  of  any  semiconductor  device 
technology.  The  RTD  model  will  be  a  semi -empirical ,  physics- 
based  model  which  will  be  scalable,  allowing  the  user  to 
accurately  simulate  RTD  circuits.  The  second  technical 
objective  of  this  project  is  to  simulate  RTD  and  RTD/FET 
circuits,  such  as  inverters,  memory  elements,  and  analog-to- 
digital  converters.  One  driving  FET  which  we  will  examine 
closely  is  the  heterodimensional  2-D  MESFET,  a  novel  side¬ 
gated  FET  which  offers  considerable  promise  for  compact,  ultra 
low  power,  multifunctional  logic  circuits.  As  part  of  this 
Phase  I  project,  we  will  also  refine  the  2-D  MESFET  model 
already  implemented  in  AIM-Spice  in  order  to  predict  nanoscale 
2-D  MESFET  device  and  circuit  operation. 
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I.  Introduction 


New  demands  for  integrated  circuits  (ICs)  are  emerging  in  areas  that  require  high  speed 
and  low  power  concurrently,  particularly  for  battery  powered  applications  such  as  wireless 
communications  and  portable  computing.  As  the  number  of  transistors  in  Very  Large  Scale 
Integration  (VLSI)  continues  to  grow  while  the  gate  dimensions  shrink,  conventional 
technology—  such  as  CMOS—  is  limited  in  terms  of  switching  speed,  power  dissipation,  and 
interconnect  capacitance.  In  addition,  cache  memory  continues  to  grow  as  more  memory  is 
integrated  locally  to  processor  cells,  increasing  the  amount  of  layout  area  devoted  to  memory 
elements.  Advances  in  semiconductor  fabrication  techniques  and  materials  processing  now 
permit  novel  nanoelectronic  architectures  with  unique  performance  capabilities  which  address 
the  limitations  of  conventional  technologies.  One  example  of  these  new  architectures  is  the 
Resonant  Tunneling  Device.  Resonant  Tunneling  Diodes,  in  particular,  feature  high  speed,  low 
capacitance,  and  material  compatibility  to  conventional  heterostructure  FETs  and  HBTs.  In 
addition,  the  NDR  of  the  RTD  load  allows  the  design  of  inverters  using  only  majority-carrier 
devices,  with  the  benefits  of  simpler  fabrication  and  high  earner  mobility  in  both  the  switching 
and  load  devices.  The  requirement  for  fast  operation  with  low-power  dissipation  in  VLSI 
circuits  suggests  that  the  NDR  of  the  RTD  is  an  ideal  inverter  load  element  since  it  exhibits  low 
current  at  digital  logic  levels  for  low  static  power  but  high  current  otherwise  for  high  speed 
switching.  Hence,  highly  functional,  compact  ultra  low  power  logic  elements  are  possible  by 
integrating  Resonant  Tunneling  Devices  with  high  speed,  low  power  compound  semiconductor 
FETs,  such  as  the  novel  heterodimensional  2-D  MESFET  [1-3]. 

A  reliable  device  model  is  a  prerequisite  for  the  design  and  production  of  any 
semiconductor  device  technology.  The  flowchart  shown  in  Fig.  1  illustrates  how  the  device 
model  is  linked  to  all  aspects  of  a  given  design  process.  For  this  DARPA  Phase  I  SBER.  project, 
we  have  implemented  a  physics-based  engineering  model  of  the  Resonant  Tunneling  Diode  in 
the  commercial  software  simulator  AIM-Spice.  We  have  also  refined  the  2-D  MESFET  model 
already  implemented  in  AIM-Spice  in  order  to  more  accurately  account  for  the  dual-gate 
geometry  of  the  device.  Together,  these  two  models  allow  for  the  simulation  of  compact,  ultra 
low  power,  highly  functional  2-D  MESFET/RTD  logic  elements. 

II.  RTD  Modeling 

There  are  two  approaches  to  modeling  Resonant  Tunneling  Diodes:  1)  numerical 
modeling,  and  2)  empirical  modeling.  To  fully  account  for  all  physical  mechanisms  associated 
with  the  optical  and  transport  properties  of  an  RTD  heterostructure,  particularly  for  high 
frequency  applications,  numerical  simulations  are  necessary.  However,  numerical  simulations 
are  computer  and  time  intensive,  so  for  applications  in  circuit  simulators,  physically-based 
analytical  solutions  are  far  more  appealing.  A  number  of  empirical  models  for  the  RTD  have 
been  reported;  however,  most  are  basically  analytical  curve  fitting  procedures  which  are  not 
based  on  the  physics  of  the  device  (see,  for  example,  [6,  7]).  For  example,  the  model  developed 
by  Robertson,  et  al,  in  [8],  which  was  used  by  Advanced  Device  Technologies,  Inc.,  in 
preliminary  simulations  of  RTD/2-D  MESFET  logic  elements,  models  the  current-voltage 
characteristics  of  the  RTD  as  a  voltage  controlled  current  source  in  parallel  with  a  voltage- 


4 


Fig.l.  Flowchart  of  integrated  circuit  technological  hierarchy  (after  [5],  p.  xiv). 


shifted  diode: 


AjV 

I  =  - - - 2  +  VBlexp{B2(V-B3)} 

\+Ay(V-V  ) 


(1) 


where  Aj,  Av  determine  the  resonant  peak  current,  is  the  peak  voltage,  Bb  B2  determine  the 
diode  current,  and  B3  represents  the  diode  voltage  shift. 
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There  are  two  notable  exceptions  to  the  analytical  approach:  the  models  developed  by 
Schulman,  et  al,  in  [9]  and  by  E.  R.  Brown,  et  al,  in  [10].  In  both  of  these  cases,  an  analytic 
SPICE  model  for  the  current-voltage  characteristics  is  developed  from  basic  principles. 
Moreover,  the  Schulman  model  is  widely  accepted  as  the  standard  model  for  RTDs.  Hence,  for 
this  DARPA  Phase  I  project  we  implemented  the  Schulman  RTD  model  into  the  circuit 
simulator  AIM-Spice. 

The  Schulman  model  is  derived  from  the  quantum  tunneling  formalism  for  the  current  in 
the  effective  mass  approximation: 


J  = 


em*feT 

27t2(/i/(27t))3 


\dET(E,V)  In 

o 


C EF-E)/(kT)  , 

1  +e 


,  (Er-E-eV)/(kT) 

1  +e 


where  the  transmission  coefficient  is  approximated  by  a  Lorentzian: 


T(E,  V)  = 


© 

1 - 1 

1 

'  2  j 

I]2*©’ 

(2) 


(3) 


Here,  E  is  the  energy  measured  up  from  the  emitter  conduction  band  edge,  Er  is  the  energy  of 
the  resonant  level  relative  to  the  bottom  of  the  well  at  it’s  center,  and  F  is  the  resonance  width. 
By  integrating  (2)  and  applying  several  reasonable  approximations  described  in  [9],  the 
following  formula  is  obtained: 


7,(V)  =  Ain 


1  +  e 


( B  —  C  +  /l]  V) 
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'7t 

2  +  atan 
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where  A,  B,  C,  D,  H,  nb  and  n2  are  fitting  parameters.  The  valley  current  contribution  is 
modeled  by  a  simple  diode. 


V 

(  n2eicT  ^ 

J2(V)  =  H\e  -1J, 


(5) 


such  that  the  total  current  is  J(  V)  =  Jx  ( V)  +  J2(  V) .  Examples  of  AIM-Spice  simulations  of  the 
Schulman  model  are  shown  in  Fig.  2. 
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Fig.2.  Examples  of  AIM-Spice  simulations  of  the  Schulman  model. 
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Fig.3.  Schematic  cross-section  of  the  Schottky/2-dimensional  electron  gas  junction 
(left)  and  perspective  sketch  of  the  2-D  MESFET  (right)  illustrating  dual  side  gates 
which  modulate  the  current  flowing  between  the  drain  and  source. 


III.  2-D  MESFET  Model 

An  accurate  analytical  model  for  the  current-voltage  characteristics  of  the  2-D  MESFET, 
based  on  the  unified  HFET  model  presented  in  [5],  was  presented  in  [1],  This  model  describes 
all  operating  regimes,  including  above  and  below  threshold,  using  a  single,  continuous  analytical 
equation.  However,  the  model  is  essentially  an  empirical  model  and  does  not  contain  geometric- 
dependencies  of  parameters  and  is  therefore  not  scalable.  As  part  of  this  DARPA  project,  we 
have  developed  a  more  advanced  2-D  MESFET  model  based  on  the  physics  of  the  device.  In 
this  new  model,  the  parameters  have  a  clear  physical  meaning  and  describe  the  geometric 
dependencies  of  the  device.  The  model  has  been  applied  to  three  gate  configurations:  (1)  the  side 
gates  are  electrically  tied  together,  (2)  the  dual  side  gates  are  biased  independently,  and  (3) 
multiple  (3  or  more)  gates  are  utilized  in  the  same  transistor.  We  have  also  developed  an 
equivalent  circuit  model  for  the  2-D  MESFET.  In  this  circuit  model,  the  device  is  split  into 
several  2-D  MESFETs,  which  are  described  using  the  unified  single  transistor  model;  hence,  the 
complete  model  is  also  valid  and  continuous  through  all  operating  regimes. 

The  2-D  MESFET  structure  is  shown  in  Fig.  3.  The  device  is  based  on  the  lateral  3-D 
Schottky  metal/2-D  electron  gas  junction.  The  Schottky  side  gates  laterally  modulate  the 
depletion  width  of  the  2-DEG  channel  and  thereby  control  the  current  between  source  and 
drain.  The  current-voltage  characteristics  of  the  2-D  MESFET  can  be  modeled  in  an  analogous 
way  to  HFETs  [5],  Hence,  the  drain  current  of  the  unified  2-D  MESFET  is  given  by  the 
expression: 


(6) 


h 


8chVds(l  +^Vds^ 


where  is  the  drain-source  voltage,  1^  is  a  unified  expression  of  the  saturation  current.  A.  is 
the  output  conductance  parameter,  m  is  a  knee-shape  fitting  parameter,  and  gch  is  the  extrinsic 
conductance  given  by: 


_ 8  chi _ 

8ch  =  1  +8chi(Rs+Rd)‘ 


(7) 


qn,\i„  .  . 

Here  gchi  =  — £ — ,  p  is  the  electron  mobility,  and  nj  is  the  electron  sheet  charge  density  per 

unit  length  given  by: 


n,  = 


(8) 


where  nmax  is  the  maximum  electron  sheet  density  per  unit  length,  y  is  a  fitting  parameter,  and 
nj  is  given  by  an  approximate  solution  of  the  unified  charge  control  model  [6]: 


n  ,  = 


2ceq*\Vth 


In 


l  +  o 


2  exp 


An 


Wth 


J  J 


(9) 


Here  r|  is  the  subthreshold  ideality  factor,  V*  is  the  thermal  voltage,  and  ceff  is  considered  as  a 
capacitance  fitting  parameter.  The  terms  Vgsl  and  Vgs2  account  for  the  dual  gate  inputs. 


The  saturation  current  is  written  as  [1]: 


7 sat 


Schi^gte^ 


(  rv*te\2\ 

1  +  8chiRs  +  /Jj  +  28chiRs  +  VL  J  ^  +  tcVgte) 


(10) 


where  £  is  the  transconductance  expansion  factor,  tc  is  the  transconductance  compression  factor, 
and  Vl=FsL  with  Fs  with  being  the  saturation  field.  The  effective  gate  voltage  swing,  Vgte,  is 
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related  to  the  gate-source  voltage,  Vgs,  and  the  threshold  voltage  Vy  as  follows: 


gte 


=  V, 


th 


V. 


1  + 


*-+  ls2  + 


r[zzzi]i 


2-i 


g* 


2r\Vth 


-1 


(ID 


where  the  parameter  8  determines  the  width  of  the  transition  between  below  and  above 
threshold.  Please  note  that  (11)  represents  the  equation  used  in  [1]  where  we  replaced  V*  with 

nv*. 


This  HFET-like  model  is  accurate  for  relatively  long-channel  devices.  In  short-channel 
devices  the  extension  of  the  effective  gate  length  caused  by  the  depletion  region  has  to  be  taken 
into  account.  We  can  model  this  effect  by  assuming  an  effective  channel  length  that  is  equal  to 
the  gate  length,  L,  increased  by  the  length  of  the  depletion  extensions,  the  maximum  value  of 

V  —V 

which  is  ddep.  We  write  ddep  =  2e  biqn~  (a  being  a  fitting  parameter)  since  this  effective 

channel  length  decreases  with  the  increasing  gate  voltage.  Therefore,  in  the  model  equations  the 
channel  conductance  should  be  replaced  by: 


_ 

8chi  ~  l  +  a ddep ' 


(12) 


We  expect  that  in  short-channel  devices  the  dependence  of  the  channel  current  on  the  gate 
voltage  will  deviate  from  linearity  because  of  the  second  term  in  the  denominator  of  the  above 
expression.  Previously,  a  different  and  empirical  expression  of  nj  was  used  to  fit  the  channel 
current  measurements  of  submicron  2-D  MESFETs  [1],  was  assumed  to  be  proportional  to 

V2gte  the  linear  regime.  No  explanation  was  given  but  it  was  actually  an  empirical  way  to 
account  for  the  dependence  of  the  channel  length  on  the  gate  voltage. 

There  is  also  a  lower  limit  for  the  effective  channel  length,  since  ddep  cannot  be  smaller 
than  zero.  The  lower  limit  of  the  channel  length  is  related  to  the  upper  limit  for  the  sheet  charge 
density,  and  this  effect  can  be  included  by  expressing  the  depletion  depth  in  terms  of  as: 


The  above  expression  tends  to  the  appropriate  limits  in  all  operating  regimes. 

The  dual-gate  2-D  MESFET  model  also  has  to  account  for  the  nonlinear  coupling 
between  the  gates  that  has  been  found  in  experimental  measurements.  This  effect  can  be 
attributed  to  two  facts:  (1)  the  asymmetry  of  the  two  depletion  extensions,  and  (2)  to  nonlinear 
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corrections  on  the  relationship  between  depletion  charge  and  gate  voltage.  This  nonlinear 
coupling  can  be  empirically  modeled  as  a  change  in  the  threshold  voltage: 

AVV=  4(VG,-VG2)2.  (14) 

In  a  dual-gate  device,  we  have  modeled  the  effect  of  the  dependence  of  channel  length  on  the 
gate  voltages  by  considering  an  average  effective  channel  length.  The  resulting  channel 
conductance  is  written  as: 


8  chi 


(15) 


,  (d,]epi  +  ^dep^) 

where  ddep  =  — • 

To  model  multiple-gate  (three  or  more  gates)  2-D  MESFETs  we  split  the  device  into 
several  dual-gate  transistors.  In  Fig.  4  we  see  the  structure  of  a  3-gate  2-D  MESFET.  For  the 
corresponding  circuit  model,  the  device  is  split  into  two  dual-gate  2-D  MESFETs,  each  one 
modeled  using  the  unified  model  described  above. 


Fig.4.  Top  view  of  a  2-D  MESFET  having  three  gates. 


We  have  validated  the  new  model  by  comparison  with  experimental  results  for  2-D 
MESFET  devices  having  different  geometries.  The  devices  were  fabricated  on  a  pseudomorphic 
heterostructure  grown  on  a  semi-insulating  GaAs  substrate.  The  parameter  extraction 
procedures  follow  basically  the  methods  reported  in  [5]  for  HFETs.  In  order  to  determine  an 
unambiguous  value  of  the  threshold  voltage,  its  extraction  is  carried  out  from  the  drain  current 
vs.  gate-source  voltage  characteristics  by  applying  the  same  voltage  at  both  gates. 
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In  Figs.  5  and  6,  we  compare  the  measured  and  modeled  characteristics  of  a  three-gate  2- 
D  MESFET  at  different  bias  conditions.  The  nominal  gate  length  is  2.0  um  and  the  maximum 
width  is  1.0  um.  In  Fig.  5  (top)  the  device  two  gates  are  electrically  tied  together.  In  Fig.  5 
(bottom),  Vgs2  =  Vgs3=  0.6  V  and  we  sweep  Vgsl,  so  the  device  is  working  as  a  dual-gate  device. 
In  Fig.  6  Vgsl  =-1.2  V,  Vgs3  =  0.6  V  and  we  sweep  Vgs2.  The  agreement  is  good  in  all  bias 
conditions,  even  near  the  threshold  voltage. 
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Fig.5.  Id  vs.  for  a  three-gate  2-D  MESFET  with  L  =  2  pm:  (top)  Vg,.  =  Vgsi  =  VgS2  =  VgS3  and 
(bottom)  Vgsl  =  Vgs3  =  0.6  V. 
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Fig.  6.  ID  vs.  Vds  for  a  three-gate  2-D  MESFET  with  L  =  2  pm  and  Vgsl  =  -1.2V  and  Vgs3  =  0.6  V. 

IV.  RTD/2-D  MESFET  Circuits 

Fig.  7  demonstrates  a  2-D  MESFET/RTD  inverter.  The  transfer  characteristics  in  the 
figure  are  for  three  different  bias  values  on  the  second  gate  of  the  device  while  the  first  gate  is 
swept  as  the  input  voltage,  demonstrating  the  dual  gate  functionality  of  the  2-D  MESFET.  (The 
voltages  at  which  the  inverter  switches  are  observed  to  vary  as  a  function  of  the  voltage  applied 
to  the  second  gate.)  Such  a  logic  element  could  be  used,  for  example,  as  a  sample-and-hold 
element  with  a  programmable  switching  point. 


Fig.  7.  RTD/2-D  MESFET  logic  element  utilizing  dual-gate  functionality  of  the  2-D  MESFET 
for  applications  such  as  a  sample-and-hold  element  with  a  programmable  switching  point. 


13 


References 


* 


[1]  W.  C.  B.  Peatman,  M.  Hurt,  H.  Park,  T.  Ytterdal,  R.  Tsai,  and  M.  Shur,  “Narrow  Channel  2- 
D  MESFET  for  Low  Power  Electronics,”  IEEE  Transactions  on  Electron  Devices,  Vol.  42, 
No.  9,  pp.  1569-1573,  Sept.  1995. 

[2]  W.  C.  B.  Peatman,  H.  Park,  and  M.  Shur,  ‘Two-dimensional  metal-semiconductor  field 
effect  transistor  for  ultra  low  power  circuit  applications,”  IEEE  Electron  Device  Letters, 
Vol.  15,  No.  7,  pp.  245-247,  July  1994. 

[3]  W.  C.  B.  Peatman,  R.  Tsai,  T.  Ytterdal,  M.  Hurt,  H.  Park,  J.  Gonzales,  and  M.  Shur,  “Sub- 
Half-Micrometer  Width  2-D  MESFET,”  IEEE  Electron  Device  Letters,  Vol.  17,  No.  2,  pp. 
40-43,  Feb.  1996. 

[4]  AIM-Spice  World  Wide  Web  page,  http://www.aimspice.com. 

[5]  K.  Lee,  M.  Shur,  T.  A.  Fjeldly,  and  T.  Ytterdal,  Semiconductor  Device  Modeling  for  VLSI, 
Prentice  Hall,  New  Jersey,  1993. 

[6]  C.  E.  Chang,  P.  M.  Asbeck,  K.-C.  Wang,  and  E.  R.  Brown,  “Analysis  of  Heterojunction 
Bipolar  Transistor/Resonant  Tunneling  Diode  Logic  for  Low-Power  and  High-Speed 
Digital  Applications,”  IEEE  Trans,  on  Elec.  Dev.,  Vol.  40,  No.  4,  pp.  685-691,  April  1993. 

[7]  Z.  Yan  and  M.  J.  Deen,  “New  RTD  Large-Signal  DC  Model  Suitable  for  PSPICE ,”  IEEE 
Trans,  on  Computer-Aided  Design  of  Int.  Circuits  and  Sys.,  Vol.  14,  No.  2,  pp.  167-172, 
Feb.  1995. 

[8]  J.  Robertson,  T.  Ytterdal,  W.  C.  B.  Peatman,  R.  Tsai,  E.  Brown,  and  M.  Shur,  “RTD/2-D 
MESFET  Logic  Element  for  Compact,  Ultra  Low-Power  Electronics,”  IEEE  Transactions 
on  Electron  Devices,  Vol.  44,  No.  7,  pp.  1033-1039,  July  1997. 

[9]  J.  N.  Schulman,  H.  J.  De  Low  Santos,  and  D.  H.  Chow,  “Physics-Based  RTD  Current- 
Voltage  Equation,”  IEEE  Electron  Device  Letters,  Vol.  17,  No.  5,  pp.  220-222,  May  1996. 

[10]  E.  R.  Brown,  O.  B.  McMahon,  L.  J.  Mahoney,  and  K.  M.  Molvar,  “SPICE  model  of  the 
RTD,”  Elec.  Letters,  Vol.  32,  no.  10,  pp.  938-940,  May  9, 1996. 


14 


Distribution  List 


Copies  Receipient 
2  Commander 

U.S.  Army  Aviation  and  Missile  Command 

ATTN:  AMSAM-RD-WS-DP-SB  (Mr.  Neal,  Tech  Monitor) 

Bldg  7804,  Room  213 

Redstone  Arsenal,  AL  35898-5248 

1  Commander 

U.S.  Army  Aviation  and  Missile  Command 

ATTN:  AMSAM-RD-CS-R 

Redstone  Arsenal,  AL  35898-5240 

1  Director 

U.S.  Army  Aviation  and  Missile  Command 

ATTN:  AMSAM-RD-WS 

Redstone  Arsenal,  AL  35898-5248 

1  Director 

Defense  Advanced  Research  Projects  Agency 
ATTN:  ETO  (Lt.  Col.  Gernot  Pomrenke) 

3701  North  Fairfax  Drive 
Arlington,  VA  22203-1714 

1  Director 

Defense  Advanced  Research  Projects  Agency 
ATTN:  OASB/ARPA  Library 
3701  North  Fairfax  Drive 
Arlington,  VA  22203-1714 

2  Defense  Technical  Information  Center 
ATTN:  Acquisitions/OCP 

8725  John  J.  Kingman  Rd.,  STE  0944 
Ft.  Belvoir ,  VA  22060-6218 

1  Director 

Defense  Advanced  Research  Projects  Agency 

ATTN:  OASB/SBIR 

3701  North  Fairfax  Drive 

Arlington,  VA  22203-1714 

1  Administrative  Contracting  Officer  * 

DCMAO  Baltimore 
ATTN:  Chesapeake 
200  Towsontown  Blvd.  West 
Towson,  MD  21204-5299 
*  transmittal  letter  only,  no  report  required  to  ACO. 


15 


